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Abstract 

We study the self-gravitating stars with a linear equation of state, P = ap, in AdS space, where a 
is a constant parameter. There exists a critical dimension, beyond which the stars are always stable 
with any central energy density; below which there exists a maximal mass configuration for a certain 
central energy density and when the central energy density continues to increase, the configuration 
becomes unstable. We find that the critical dimension depends on the parameter a, it runs from 
d — 1 1 .f429 to f0.f29f as a varies from a = to f . The lowest integer dimension for a dynamically 
stable self-gravitating configuration should be d = 12 for any a € [0, 1] rather than d = 11, the latter 
is the case of self-gravitating radiation configurations in AdS space. 
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I. INTRODUCTION 

Self-gravitating configuration is a subject of long-standing interest in general relativity since its thermody- 
namic behavior is quite different from that of a usual thermal system without gravity, due to the attractive 
long-range, unshielded nature of gravitational potential. For example, it is well known that the canonical 
ensemble is not defined in asymptotically flat space. This is because having thermal radiation at constant 
temperature at infinity is not compatible with asymptotic flatness. One can avoid this problem by enclosing 
the system in a box, which is unphysical, or by working in anti-dc Sitter (AdS) space which needs not any un- 
physical perfectly reflecting walls at finite radius, the rising gravitational potential in AdS space, plus natural 
boundary conditions at infinity, acts to confine whatever is inside. On the other hand, due to the conjecture 
of AdS/CFT correspondence which says that string theory/M theory on an AdS space (times a compact 
manifold) is dual to a strong coupling conformal field theory (CFT) residing on the boundary of the AdS 
space, over the past decade, a lot of attention has been focused on AdS space and relevant physics. Further, 
Witten [2j argued that thermodynamics of black holes in AdS space can be identified to that of dual strong 
coupling CFTs. Therefore one could study thermodynamics and phase structure of strong coupling CFTs by 
investigating thermodynamics and phase structure of AdS black holes. It is well-known that thermodynamics 
of AdS black holes is quite different from that of their counterparts in asymptotically flat space and that there 
is a Hawking-Page phase transition between large black hole and thermal radiation in AdS space [3( . Therefore 
it is also of great interest to study self-gravitating configurations in AdS space and their thermodynamics. 

Self-gravitating radiation gas has been investigated thoroughly. Sorkin, Wald and Zhang 4[ have studied 
the equilibrium configurations of self-gravitating radiation in a spherical box of radius R in asymptotically 
fiat space. It was found that for locally stable configuration, the total gravitational mass of radiation obeys 
the inequality M < ^ max -R, where /i max = 0.246. In AdS space, Page and Philips [5| examined the self- 
gravitating configuration of radiation in four dimensional space-time. The configuration can be labeled as its 
mass, entropy and temperature versus the central energy density. They found that there exist locally stable 
radiation configurations all the way up to a maximum red-shifted temperature, above which there are no 
solutions; there is also a maximum mass and maximum entropy configuration occurring at a higher central 
density than the maximal temperature configuration. Beyond their peaks the temperature, mass and entropy 
undergo an infinite series of damped oscillations, which indicates the configurations in this regime are unstable. 
The self-gravitating radiation in five dimensional AdS space has been studied in [6j (see also 7|) with similar 
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conclusions. Recently, Vaganov 

8] 

and Hammersley [9j independently discussed the self-gravitating radiation 
configurations in higher dimensional AdS spaces. They found that in the case of 4 < d < 10, the situation 
is qualitatively similar to the case in four dimensions, while d > 11, the oscillation behavior disappears. 
Namely, there is a critical dimension, d c , a ds = 11 (very close, but not exact), beyond which, the temperature, 
mass and entropy of the self-gravitating configuration are monotonic functions of the central energy density, 
asymptoting to their maxima as the central density goes to infinity. The equilibrium configurations of self- 
gravitating radiation gas in AdS space are quite different from those of their counterparts in asymptotically 
flat space. 

Recently Chavanis investigated spherically symmetric equilibrium configurations of relativistic stars with 
a linear equation of state P — ap in flat space [10j. The equation of state P — ap says the pressure P 
is proportional to the energy density p, where a is a ratio coefficient. Theoretically, different coefficients 
correspond to different stars, although it is impossible to define names for all the stars corresponding to each 
coefficient. The speed of sound is given by (dP/dp) 2 c = a 2 c, where c is the speed of light. Thus causality 
requires a 2 < 1. As a result, < a < 1 is considered only. One could identify relativistic stars by an equation 
of state with their different ratio coefficients. In the Newtonian limit a — ► 0, it returns the classical isothermal 
equation of state. On the other hand, the case with a = l/(d — 1), where d is the dimension of space-time, 
corresponds to a gas of self-gravitating radiation; the core of neutron stars or so called "photon stars" where 
the pressure is entirely due to radiation. A gas of baryons interacting through a vector meson field for the 
case of a = 1 is called the "stiffest" star. Chavanis found that the structure of the system is highly dependent 
on the dimensionality of space-time and that the oscillations in the mass-central density profile disappear 
above a critical dimension d c .fi a t(a) depending on the coefficient a. Above this dimension, the equilibrium 
configurations are stable for any central density, contrary to the case d < d c ,fi a t- For Newtonian isothermal 
stars (a — * 0), the critical dimension d c ./;at(0) = 11. For the stiffest stars (a = 1), d c ./iai(l) = 10 and for a 
self-gravitating radiation(a = l/(d— 1)), d c .fi a t — 1 + 9.96404372... very close to 11. The oscillations exist for 
any a G [0, 1] when d < 10 and they cease to exist for any a G [0, 1] when d > 11. 



In a previous paper 11(, we have studied self-gravitating radiation configurations with plane symmetry in 
AdS space and found that the situation is quite different from the one in spherically symmetric AdS space. 
In the present paper we continue this study and consider the self-gravitating configurations with a linear 
equation of state P = ap in higher (e? > 4)-dimensional spherically symmetric AdS space, in order to see the 
dependence of the critical dimension on the parameter a and to see whether there is any essential difference 



between the configurations in flat space and the configurations in AdS space. 

The organization of the paper is as follows. In the next section, we give a general formulism to describe 
the relativistic stars with a linear equation of state P — ap in AdS space. The numerical results are given in 
Sec. III. The Sec. IV is devoted to the conclusions. 

II. RELATIVISTIC STARS WITH LINEAR EQUATION OF STATE IN ADS SPACE 

Consider a d-dimensional asymptotically AdS space with metric 

ds 2 = -e 2S{r) h(r)dt 2 + h^irjdr 2 + r 2 j ij dx i dx j ', (1) 

where S and h are two functions of the radial coordinate r, and 7y is the metric of a (d — 2)-dimcnsional 
Einstein manifold with constant scalar curvature (d — 2)(d — 3). We take the gauge linv^oo <5(r) = 0, and 
rewrite the metric function h(r) as 

, , , , r 2 l6irGm(r) 

h ^ = l+ p- ( d -mr^ (2) 

where I denotes the radius of the AdS space with cosmological constant A = — (d — l)(d — 2)/2/ 2 , is the 
volume of the Einstein manifold, and m(r) is the mass function of the solution. In this gauge, the total 
gravitational mass of the solution is just 

M = lim m(r). (3) 

r — >oo 

The Einstein field equations with the cosmological constant and energy-momentum tensor T^ v are 



1 (d-l)(d-2) 
n^v - -jg^vti ^ 9^ = (4) 



The metric function S(r) and the mass function m(r) satisfy the following equations 

S'( r ) - 8nGr (T* - T r ) 

m'(r) = -Vtr d - 2 T\. (5) 

The stress-energy tensor of perfect fluid is = (p + P^U^XJu + Pg^, where is the four- velocity of the 
fluid. For self-gravitating radiation, the equation of state obeys P = p/(d — 1). For linear relativistic stars, 
P = ap with a € [0, 1]. Thus the equations in ((SJ reduce to 

m'(r) = VLr d - 2 p. (7) 
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FIG. 1: The self-gravitating radiation in AdS space: The mass of the self-gravitating radiation configuration versus 
the central energy density from dimension d = 4 (top curve) to d — 12 (bottom curve). 

In this paper we are particularly interested in the relation between the total mass of configuration and the 
central energy density. To integrate ((6|) and (|7|), we make a scaling transformation as follows, 



r— >lr, p—*l p, m(r)—*l m(r) 



(8) 



so that r, p and m become dimensionless. In the numerical integration, we will adopt the units 8ttG = 1 and 
1 = 1, and rescale the mass function as 

16nGm(r)/(d - 2)fi -> m(r). 

In that case, the gravitational "mass" M in the plots in the next section in fact is the gravitational mass 
density, lQirGM/{d — 2)fl, of corresponding relativistic star configurations. From the conservation of the 
energy-momentum tensor, one can derive the energy density p(r) satisfies the following equation 

dp (a + \)pr ( d — 3 d — 1 d — 3 2ap 



dr 2a V t 2 h ( r ) r 2 /i(r) (d - 2)h{r) J ' ^ 

Further, note that the equations ([7]) and ([9]) are singular at r — 0. To avoid this, in the numerical calculations, 
we will start the integration from r = e = 1CP 5 to r = L = 100. Obviously, the accuracy of the numerical 
calculations depends on the values of e and L. 



III. NUMERICAL RESULTS 



To be more clear to discuss all kinds of critical dimensions of relativistic stars with a linear equation of state 



in AdS space, let us first revisit the self-gravitating radiation case. In 



Hammersley found that in the case 



of 4 < d < 10, there exist locally stable radiation configurations all the way up to a maximum mass, beyond 



10 8 ~ 2d (d-3) 2 M 
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FIG. 2: The case of a = 0.5 in AdS space: The mass of relativistic star configurations versus the central energy density 
from dimension d — 4 (bottom curve) to d — 12 (top curve). 

their peaks the curves undergo an infinite series of damped oscillations, which indicates the configurations 
in this regime are unstable. However, with d > 11, the oscillation behavior disappears. The configurations 
turn to be monotonic functions of the central energy density, asymptoting to their maxima as the central 
density goes to infinity. Namely, there is a critical dimension. By numerical analysis, Hammersley [9] found 
a semi-empirical model 

log Pe ^0.50d+ ' -2.20, (10) 
Vll.O — a 

which gives a critical dimension d c a( i s — 11, and a similar conclusion was also reached by Vaganov 8| 
independently. We reproduce the result in Fig. [1] 

Now we turn to the case of relativistic stars with a linear equation of state P = ap. Due to the requirement 
of causality, we consider the cases with < a < 1 . In Fig. [51 we plot the mass of relativistic stars configurations 
versus the central energy density from dimension d = 4 to d = 12 for the case of a = 0.5. Note that in Fig. [5] 
we rescale the mass M with scale 10 8 ~ 2d (d— 3) 2 in each dimension for better displaying all curves in one figure. 
We notice that for a certain value of a, the configurations of different dimensions from 4 to 12 are similar to 
the case of self-gravitating radiation stars. For the case of lower dimensions, there exist obvious oscillations 
which can be seen from Fig [3l As for the higher dimensions, the oscillations become weaker and weaker, 
and finally beyond some critical dimensions, configurations become monotonic functions of the central energy 
density. Fig. [4] is the case of d = 11. Note that in Fig. [3] and Fig. [4] each curve corresponds to each relativistic 
star with coefficient a which undergoes a continuous variation from 0.01 to 1.0 with step 0.1, and also we 
rescale the mass density M to M/M max in each case of coefficient a for better comparing those different cases. 

As a matter of fact, it is not really accurate to determine the critical dimensions only by the naked eyes to 
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FIG. 3: The case of d=4 in AdS space: The mass of relativistic star configurations versus the central energy density 
from the coefficient a = 1 (left curve) to a = 0.01 (right curve). 

M 




FIG. 4: The case of d=ll in AdS space: The mass of relativistic star configurations versus the central energy density 
from the coefficient a = 1 (left curve) to a = 0.01 (right curve). 

watch these curves. To be more accurate to determine the critical dimensions, we take the numerical analysis 
approach. In Fig.[5]one can see clearly that for each curve there exists a saturation point, log p c (the red dot in 
the figure), which is the location of the first local maximum. The saturation point moves towards the right side 
as the dimension d increases. Beyond the saturation points, the curves undergo an infinite series of damped 
oscillations. When the dimension increases, the oscillations become weaker and weaker, while the saturation 
point log p c becomes larger and larger, and finally, beyond a critical dimension d c . a ds, the saturation point 
log p c goes to infinity. So we can determine the critical dimension by analyzing the variation of the saturation 
points. 

Some data of saturation points log p c are listed in Table Q] Obviously, log p c depends on the coefficient a 
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FIG. 5: The self-gravitating radiation in AdS space: The mass of the self-gravitating radiation configuration versus 
the central energy density from dimension d — 4 (top curve) to d = 10 (bottom curve). The red dots are the saturation 
points log p c . 

and dimension d, so namely logp c (a,d). We use a formula like (|10j) and obtain the critical dimensions d c (a) 
for different coefficient a, which are listed in Table [Hi As can be seen from Table[l] in the case of d = 11 some 
data singularity appears. So we only use data from d = 4 to d = 10 to fit. 



TABLE I: Saturation Points log p c (a,d) 



a 


d 


=4 


d 


=5 


d=6 


d=7 


d 


=8 


d=9 


d=10 


d=ll 


0.01 


3 


90 


3 


60 


3.90 


4.20 


5 


10 


6.00 


8.10 


21.00 


0.1 


3 


00 


3 


00 


3.30 


3.90 


4 


80 


6.00 


8.40 


21.90 


0.2 


2 


40 


2 


40 


3.00 


3.90 


4 


80 


6.30 


8.70 


oo 


0.3 


1 


80 


2 


10 


3.00 


3.90 


4 


80 


6.30 


9.00 


CO 


0.4 


1 


50 


2 


10 


3.00 


3.90 


5 


10 


6.60 


9.60 


oo 


0.5 


1 


20 


2 


10 


3.00 


3.90 


5 


10 


6.60 


10.20 


oo 


0.6 


1 


20 


2 


10 


3.00 


4.20 


5 


10 


6.90 


11.10 


oo 


0.7 


1 


20 


2 


10 


3.00 


4.20 


5 


40 


7.20 


12.30 


oo 


0.8 


1 


20 


2 


10 


3.30 


4.20 


5 


40 


7.50 


14.10 


oo 


0.9 


1 


20 


2 


10 


3.30 


4.20 


5 


70 


7.50 


15.60 


oo 


1.0 


1 


20 


2 


10 


3.30 


4.20 


5 


70 


7.80 


16.50 


oo 



TABLE II: Critical Dimensions d c . a ds{a) 



a 


0.01 0.1 


0.2 


0.3 


0.4 


0.5 


d c 


ads (a) 


11.1061 11.0952 


10.9989 


10.9254 


10.7465 


10.6333 


a 


0.6 


0.7 


0.8 


0.9 


1.0 




ads(a) 


10.4928 


10.3707 


10.2602 


10.2023 


10.1763 



TABLE III: Numerical Fitting Results 



deads (a) 


linear 


second-order 


third-order 


average 


low bound 


10.1063 


10.1048 


10.1763 


10.1291 


up bound 


11.1601 


11.1591 


11.1097 


11.1429 
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11.2 



11 



10.8 
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FIG. 6: The critical dimensions d c . a d s (a) versus the coefficient a. The blue hollow dots are the scattered data of the 
critical dimensions d c , a d s (a). The red dashed curve is the case of linear fitting. The blue dashed curve is the case of 
the third-order fitting. The black thick curve is the case of average fitting. 

In Fig. [5] we plot the scattered diagram of the critical dimensions d c , a ds{o-) versus the coefficient a. It can be 
seen that the value of d c . a d s (a) becomes smaller and smaller as the value of the coefficient a goes from to 1. 
So d c a ds{o) can be considered to be a monotonic function of the coefficient a. The red dashed line in Fig.[S]is 
the result of a linear numerical fitting for these scattered data. It could be seen that for all kinds of relativistic 
stars with a linear equation of state P = ap (a £ [0, 1]), the critical dimensions d c . a d s (a) £ [10.1063, 11.1601]. 
Obviously, the lower and upper bounds of d c . a d s {o) for this kind of fitting are a little wider because these 
scattered data are not really linear related. To obtain more accurate results, we need to adopt non-linear 
fitting which could better fit the scattered data. The results of the second and third-order fitting are listed 
in Table IIII1 The blue dashed curve in Fig. [5] is the case of the third-order fitting. Note that the curve of 
the second-order fitting is not plotted in Fig. [6] because it is overlapped with the most of the linear fitting 
curve. It can be seen from Fig. [5] that the curve of the third-order fitting can better fit the scattered data. 
Namely, in AdS space for all kinds of relativistic stars with a linear equation of state P = ap (a £ [0,1]) 
the critical dimensions d c . a d s {a) £ [10.1763,11.1097]. The minimum for the stiffest star in AdS space is 
rfc.ads(l) = 10.1763 and the maximum for the Newtonian isothermal star is d c . a ds(0) = 11.1097. In Fig. [7] 
we plot the saturation points log p c versus the dimension d for each relativistic star (corresponding to each 
coefficient a ) from a — 0.1 (right blue curve) too= 1.0 (left red curve) in the range 9.4 < d < 11.4. The low 
and up bounds of the critical dimensions d c . a ds{o) £ [10.1763, 11.1097] are also plotted in this figure (the blue 
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FIG. 7: The saturation points p c versus the dimension d from a = 0.1 (right blue curve) to a — 1.0 (left red curve). The 
blue dashed curves are the low and up bounds of the critical dimensions d c . a ds(a) £ [10.1763, 11.1097], respectively. The 
black thick curves are the low and up bounds of the critical dimensions d c . a ds.aver(p) £ [10.1291, 11.1429], respectively. 

dashed curves). One can find that for the stiffest stars (a — 1), the left red curve increases steeply and finally 
up to a critical dimension d = 10.1763. However, for Newtonian isothermal stars a — > 0, the right blue curve 
changes more slowly, up to the critical dimension d = 11.1097. 

By an analysis of asymptotic behavior, Chavanis 10] found that the critical dimension runs from d = 10 



to 11 when the coefficient a decreases from a = 1 to in flat space. Namely the difference of dimension 
Ad/jnt = d c .fiat(0) — d c .fiat(l) = 1. In our case, namely in AdS space, the critical dimension runs from 
d = 11.1097 to 10.1763 with Ad a d s = ^c.ods(O) — rfc.ads(l) = 0.9334 when a varies from to 1. If taking the 
averaged values, then one has Ad ads = d c , a d s . a ver(0) — d c , a ds.aver{^) = 11-1429 — 10.1291 = 1.0138. There 
exists some difference between the flat and AdS cases. Therefore we guess that for all kinds of relativistic 
stars with a linear equation of state regardless of in AdS and fiat spaces, the critical dimensions d c . a ds(a) 
and d c ji at (a) vary just within 1 dimension when a changes from to 1, and the inaccuracy between the two 
cases is very small, within 1.38%. Therefor the critical dimension for a stable relativistic stars with a linear 
equation of state is 12, rather than 11, which is the case of self-gravitating radiation configurations. 



IV. CONCLUSION 



There exists a critical dimension for self-gravitating configurations in general relativity, beyond which the 
configurations with any central energy density are always stable; below which there exists a maximal mass 
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configuration for a certain central energy density, when the central energy density increases, the configuration 
becomes unstable. In this paper we studied the self-gravitating configurations (relativistic stars) with a linear 
equation of state P = ap in AdS space, where a is a constant parameter within a £ [0, 1]. We found that the 
critical dimension depends on the parameter a, it runs from d = 11.1097 to 10.1763 in the third order fitting 
as a varies from a = to 1. The result is a little different from the case in flat space. It runs from d — 11.1429 
to 10.1291 if one takes the averaged fitting. In that case, it is very close to the case in flat space within the 
inaccuracy 1.38%. Therefore it is of interest to study the self-gravitating configurations in de Sitter space and 
to see whether there exist any differences among the three cases. 

In 8j and 9j , it was found that the critical dimension of self-gravitating radiation configurations is d = 1 1 . 



Combing the results obtained in ll| and in the present paper, we see that in fact, the dimension d = 11 
might not have any other special physical meaning, except for the stability divide. Note that the fact that no 
special happens for thermodynamics of AdS Schwarzschild black holes in d > 4, it is of interest to understand 
the meaning of the existence of the critical dimension in the AdS/CFT correspondence. 
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